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2) ‘..'hich avulUs 


icr. of cluster calculPtions using yuih tlie 
p K*-— < '•'•’"''vtftionc.l rccr.r.iq«^ (Jwuijj" 
s i,.uLj.'i.u- uiu appi'oKir.acion, we nave 
investigneed the energy levels in clusters of nickel atoras and 
compared the results. The point group symmetry of the cluster 
has significant effect on the energy of levels nominally not 
occupied. This will influence the electron transfer process 
during chemisorption. The SSO technique permits us to treat, as 
a cluster, the approaching atom or molecule plus a small number 'of 
nickel atoms. Specifically, we have calculated that molecular 
levels become more negative in the 0 atom, as well as in a CC 
molecule, as the metal atoms are approached. Thus, electron trans* 
for from the nickel and bond formation will be facilitated. This 
new type of result is of importance in understanding chemisorp- 
tion and catalytic processes. 


^Research supported by a grant from the National Aeronautics and Space Adminis- 
tration and a contract with the U, S. Atomic Energy Commission, 
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The subject: oC chentisorption has long been a subject of scientific interest- 

since the distinction bctv;ocn it and physical adsorption v,’as made before 1930. 

* * * 

Since this paper vill be ^>5VOtcd to calculation of a chemical physical nature, it 

is v/oll to cite so:.to oc the criteria by v?hich chemisorption is recognized^^ 

o\ 

(Rioting from Trapncll*'' , '’a);. Adsorption is localized and tahes place through 
S 

collision of gas molecules with vacaat sites of the substrate, b), each site can 

♦ 

accomodatd one and only one adsorbed particle molecule and c). the energy of an 
adsorbed particle is [substantially] the same anywhere on the surface, and is 
independent of the presence or absence of nearby adsorbed molecules.*' 

The experimental developments have been very 'significant and the migration 

31 

of species has been demonstrated by Ehrlich , a number measurements on the 

4N 5) 

kinetics of activated adsorption have been made ' and Gomer and his Colleagues 
have observed surface structures directly. The disruption of su'^rface layers as 
observed by Germer and MacRac^^ for the chemisorption of hydrogen on nickel is 
particularly significant since it stimulated many studies which showed that 
reconstitution apparently occurs in many surface systems. 

In contrast, the theoretical analysis has not kept pace in this author's 
opinion. 

Various hypotheses have been advanced to explain the relatively strong bond 
which Is formed betx^cen the adsorbate and the substrate. These have perforce 
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' Uc&cly dealc with phanoir.cno Logical aspocca of tho adsorption process. It is 
vorthvhilc to sununarizc the various approaches token before turning to the 
present method. calculated adsorption bond energies by taking the 

difference of the bund energy of (for nccaL-hydrogcn bond) and (the 

t 

dissociation energy of a molecule.) He also included the difference in the 

O \ 

electronegativities •/, proposed by Pauling^'. Eley included a dipole energy, 
again obtained from Kvan v;ith a variety of approKimations and improvements 

the energies arc not in good agreement with e:;pcrimcnt. The difficult point is 

. * • 

to adequately estimate die' term i.'hcn only s. few ’’bonds" arc broken. 


9 ) 

Trajjnell ' has obser/cd that the heats of adsorption for a variety of gases 
such as O^, ^2, CO, CO^, etc. generally ciacreasa for the series of transition 
metals as .one goes from Ti to h'i and Cu. This trend corresponds to a decrease 
in the number of available unoccupied d-levels. So, the d electrons participate 
in sore 'Significant way and many of the tlioureticai approaches to be discussed 
below nave their greatest difficulties in dealing v/ith such electrons. Further, . 

the misleading statements exist in the literature which suggests that these 

«» 

electrons are very localised in transition metals -and they cannot' participate in 
bonding. Waber^^^ has presented the redial distribution of d-li’xe crystal 
orbitals (obtained by the Augmented Plane Wave Method) and it is clear that the 
d-electrons have appreciable probability of being found at the surface of a 
Wigner-Seitz sphere. . . 


Pauling^^^ has ‘ introduced the idea of resonating covalent bonds (such as dsp) 
and this has brougtit the thinking of many scientists--dcspitc the fact that they 
may not agree with all the details) to the point that the d-character of possible 
bonds are included in their thinking and explanations. 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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Another type of treatment assumes considerable electron .transfer and the 

* 

4 > m 

formation of M A bonds, then the problem of chemisorption fs reduced to making; 

calculations of the fiorn-Mnycr type. While such an approach has the merits of 

simplicity, it tends to oversimplify, to ignore Image potentials (due to charge 

accumulation on the surface of a continuum.) In these attempts, the metal la 

12 ) 

not carefully included in the system, Dennett and Falicov attempt a more 
careful analysis but in the end, resorted to treating tungsten as a free-electron 

t * 

metal. , 


A variety of quantum-mechanical approaches have been used. These started 
with and v;ere more carefully studied by Higuchi and Eyring^^\ They 

treat the surface compic:-: as a diatomic molecule. The diagonal element of the 
Hamiltonian (v/hich v;ill be discussed below in terms of the present method) 


they write 


“ V^e 


ORIGINAL PAGE IS 
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where A^ Is the electron affinity and I is the ionization potential. Here Is 
the radius giving the minimum. Again, they resort to substantial approximations 
and the interpretation in terms of partially ionic character of the bond formed 
in the complex, does not significantly modify the valuable ad hoc approach taken 
by Eley and Pauling, 


Another popular type of treatment is to use the Valence Bond approach 

particularly in terms of the angular parts of atomic orbitals. Note that the 

potential field of the metal docs not give rise to spherically synmetrlc 

potentials. Thus, the degeneracy of certain atomic levels is lifted as it is in 

crystal field theory ."ind Ugand field theory. Classical treatments of that type 
X 5 ^ 1 Q 

. arc by Coodenough*'^^ and Orgcl \ These authors deal with primarily atoms 
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<or Ions) surrounded b/ a uniform or symmotrlcol Arrangement of Ions, Belkhausan * 

treats the more complex problem of electron clouds surroundlojo' "He ligands 

rather than point charges. K more detailed attempt to take > account is 

the lou^v symmetry of the arrangement of substrate atoms (on various planes) 

which was recently published by Rao and Waber^^^. Their specific assumption was 

that as far as an approaching atom was concerned, it experienced a non-uniform 

potential from the nearby substrate atoms as though they had a net unsatisfied 

charge in the vicinity of 0.1 of an electron. ORIGIt^AI. PAGE IS 

OF POOR QUALITY 


To this collection of quasi-chemical approaches, the Important review 

ig\ 

by Schrikifer , should be added, to put the more physical aspects into proper 

20 '’l^ 

perspective. This contains discussions to the work’of Grimley , Weycicchowski" :, 
22 ) ?. 3 ; 24 ) 

Gadzuk , Kc^tvs*’ and Plummer and Rhodin . Unfortunately the excellent recent 

*>5 26 ) 

work by Lange and Kohn" ’ is not refcircnced. Space limitations cause us to 

ou MU ^ xHtcrested reader xs exreccec to 

these primary sources to amplify his understanding of recent quantum mechanical 
approaches. Schriefer's review puts the various treatments into a framework 

m 

where they can be judged against one another. 


Several recent attempts have been made to clucid^'^te the mechanism in quanti- 
tative form using molecular orbital approaches. Of this genre, the work by 
Messmer ^ ^ on the adsorption of various gases on layers of graphite could 

be cited aa an example of one molecular orbital approach to surfaces. To avoid 
the calculation of certain integrals, he resorted to an approximation known as 
Complete Neglect of Differential Overlap (CNDO.) Another impurcant recent develop- 
ment was the adaptation of the multiple scattering Xa (MSXa) method of molecular 
calculation dcv<’llopcd by Rclth Johnson, John Slater and their collaborators^^*' 

For clusters, this method characteristically relics on treating the potential 
around cad*, atom site in a special, way -- tl'.at is, inside each sphere of radius 
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of Che ppccnclaL ii subscanclally acomlc (cencraL field epproxitnacion,) 

However^ tn the region oucsidc r^ the poccnelaL is spherically averaged oc a * 

conscanc value. This can be described as a "nuffin-cin" pocential since Che 

radial regions cencerod at each of che and nuclei* in a clusCer may be 

in concaec buc do not overlav'i i*e. , r^ rj £ d^^i where d^j is Che inCeraComic 

distance. The values of and Vj are generally chosen Co be similar Co Che 

3n 32) 

atomic (or ionic) radii given cither by Slater , Waber and Cromer or 

Pauling^^\ A very recent development by Johnson independently by 

’ * 55 ) r' 

Liberman ^ is to pennit the spherically averaged cKtcrnal regions to overlap 

so Chat d_. v/ill be less than the sum of r and r.. Danese and Conklin^^^ 

IJ 1 M 


have bepn investigatlns this method. 


ORIQINAI. PAGE 
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The present method rotnins most of the very desirable features of the (MSXa') 
method of molecular calculation but further relates this condition of spherically 
smoocnea ootentials. It has evolv*^<< frnm » nnober of •'♦•.iidi''®' ""rrio'^ 

out by Ellis and Painter^^^ using the Dis rete Variational Method (DVM) to calcu- 
late observable quantities of molecular clusters. 


In DVM one presumes that they can obtain a very’ good representation of the 
values of the collection of Important interatomic integrals by evaluating their 

37 3 

integrands at a discrete set of points. For this purpose^ Ellis and his coworkers ’ 

used a basis set of Caussians and Slater>typc orbitals in Che regions outside 

r. and r, but solved the atomic Hamiltonian inside the individual spherical 
I J 

regions surrounding the 1~" and J— ion cores. To emphasize this difference in 
Che two approaches, one notes that Johnson and coworkers used spherical waves 

/TAN 

(similar to those used in the Korringa-Kohn-Rostackcr'’'’''' (KKR) method for band 

40) 

structure calculations) in the cxtrasphcrical (or interatomic) region. Painter 
Investigated* the use of a linear, co'nibinatiun of KKR orbitals in connection with 
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OmOINM. * 

the Discrete Variational Method for metallic bonds. OF POOR QUALITY 

<> 

At essentially the same time, Averlll and Ellls^^^^ and Painter and M. Stock^^^ 
recently proposed usins ^ linear combination of atomic orbitals. Because the 
orbitals for excited states tend to be diffuse, they Imposed a high potential 
barrier at a reasonably large distance so as to localiec the Rydberg orbitals 
but not disturb the minimal basis set of hound atomic orbitals. For this reason 
they labelled their method Single Site Orbital <SS0). Because of Close collabora- 
tion over the years, these tv/o treatments arc very similar. 

It ^should be pointed out that: the (MSX:^) method the cluster of atoms is 

confined by a larger sphere. Thus it docs not readily lend Itself to treatment 
« 

of on atom or molecule approaching a cluster. In contrast, restriction is lifted 
since the SSO was specifically designed to treat the angular variaeion of the mol- 
fewultt*. li’i tliu bcuv;een and outside the spheres, wniie keeping 

the excited orbitc.ls In the basis set "well behaved." 


« 

Mothod og CompueoClon 
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The basic pcoblcsn in all such atomic and molecular structure problems is to 
obtain clgon-values and qii'onfunctions for a one-electron Hamiltonian 3C where 
vte') represents some approximation to the potential an electron experiences in 
the molecule. Such an eisenfunction may be written as 


Xi (*) ■ Sxj cs> c,, 


( 1 ) 


where th|p functions Xi (r) arc symmetrised linear combinations of single site 
orbitals % 


Xj <5) " 5,„ “aj 


( 2 ) 


where the the radial and angular functions centered at the 

u~ ion core and n, i amd n are the appropriate quantum numbers. There is not 

sufficient space to show how the symmetrization coefficients can be obtained by 

group theory--it is sufficient to indicate that j identifies a particular basis 

function for the irreducible representation. The radius r^ is the distance 
CK 

'between the nucleus and the radius vector^. 


The matrix elements which need to be calculated are 


-<XiMxi> 


(3) 


S|J -<xt|xj> 


(4) 


3 
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Th* iliicinsuialilnc CoaCucc of tho DUcrcM Vorloelonal Method in contrast 
to the RaylcigU-hitz method, is that such matrix olomcnts arc evaluated as 


j<ii X, [-7^ + >1, 


where the three components (in real space) are detei mined randomly. 

Hero is an appropriate weighting factor applied at each of the K points. 


The spacial density of points need not be uniform. In facti Averill and 
Ellis usje a distribution which is similar to the two parameter Fermi distri- 
bution of protons in a nucleus, that is • 


■ 


w 


4 Tt Vy -{l + exp !■ j. 


< 6 ) 


where is a nonnallzation constant. The probability of a point occuring in a 

spherical shell of radius r Is A /2 when r . The distribution is essentially 

u u u u 

uniform at the value A^ until is close to and then approaches zero rapidly 
with a "decay” constant or "e-folding" distance of 

Since these distributions of sample points are centered about each ion core, 

* * . 

Averill and Ellis use a linear combination of these 

»('> • 5 “uCji 

t 

where t^ is an adjustable sampling factor for the u~ atom. They chose 

B^*>1 and R* equal to the atomic radius. The total, number N of points at viiich 

>1 

» 

y 


thtt Incogrands aro sampled influence cho accuracy of finding cigenvralufis. 


To obtain the Single Set Orbitals, they chose 


V '^(r) 

fso' 


v'aCO-V^ 

0 


r 


r >R 


( 8 ) 


and solved the equation 




V 
Attin 




0 ) 


for the atomic like eigenfunctions. 


iTt,« 5*f »jound atomic 


minimal, but wuO tw lucluuu 


interatomic interactions. However, reasonable Rydberg states can be found by a 
Judicious choice of the well depth and the well radius without affecting 
the bound orbitals appreciably. 


' ORIGINAL PAGE IS 
OF POOR QUALITY 
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Th« rasuLcs cbealned to-4at;tt ate of two types: (a) the oner^ levels in small 
clusters of nickel atoms were found by a variety of different methods and for two 
spotial conformations and T^j— these results are compared below— and (b) the 
shift in the electronic levels of an oxygen atom and of a CO noleeule as it 
approaches such u cluster of ^li atoms* the latter result is of significance for 
chemisorption and for thc^. react ions between an atom and an intorfacei 


In fis, 1, the eisenvalucs obtained after the first iteration are compared 

for tw’o *dif|crent arrangements of five nickel atoms. The free atom values are for 

8 Z 

3d and 4s, assuming an initial configuration of 3d 4s , Because of the higher 


symmetry in the T 
levels than for D 


cluster, thci <2 are fewer (but more degenerate) molecular 
this is readily seen for the hybrid excited levels, the 

1% If AM i*U#*M^* * 1 ♦ Iia S 

iir ww*»w owW<swatW mm msowssss *sa # *»09 


A very important variation of the valence eigenvalues with method of calcula- 
tlon is shotm in Fig. 3. In the first group, are the results for the MSXa 
’•lluscatel" progrmn developed by Johnson al. The second group is for a 
nodlficatlon of the present method (DVM-SSO) in which spherical averaging of the 
interatomic potential has been Imposed to facilitate comparison with the MSXaf 
results, the effect of relaxing this arbitrary condition in SSO calculations is also 
shown, the shifts in the levels near *0.7 Ry are interesting, these levels were 
obtained with N«S0O sampling points. In the first three panels of Fig. 3b, the 
molecular levels obtained for 1000, 2000 and 4000 sampling points in the 
Dlophantinu intc;',rat:lon scheme of DVH. ,e is surprising that the principal 
shifts occur in BIG level at the top of the drawing as well as for the four 

t 

hydridised levels below >0.7 Ry. An addicional set of results ’obtained for an 
Increased intecacomLc separation in the conformation are show in the last 


II 


p«nol of Fic. 3b, T!ic low lying l«vcU were not plotted for volenco tlectronl 


The graphff nhowin);;' »lniUr offoecs on levels derived from 3p end 3s states 
have been plotted but they ere not Included. 

Since no efforts heva been made to carry out the DVlt-SSO results towards 
solf-consistoney, it is important to investigate the effect of the successive 
iteration on the Huscatel oiscnvalues* This is presented in Fig. 4. While some 
shifts do occur with the approach towards self-consistency* the excited levels 
Item insensitive. The major changes occur in the ordering of the strongly 
intcracblng Levels derived from the 3d atomic state. Large shifts are also appar 
ent* But after 15 cycles, the levels do not vary greacly from those obtained 
after the first iteration. This data is of interest tc both a practitioner 
in the field and to cKperimontalists wishing to interpret spectra, but has little 
influence on tritr* general lo *''y ^ir'*t-p?.cs 

Td and calculations presented in earlier figures. 

These results on small clusters may be compared with the results soon to 

43 ^ 

be reported by Diamond, Messmer, Knudson and Johnson whom studied clusters 
containing 13 or more atoms of nickel and copper. There, more states will be 
found because of the difference in coordination numbers for atoms in the 
suboctahcdral and icosEhedral clusters. . 

< 

The two portions of Fig. S arc composite drawings. On the left -hand 

I 

side, the density of curve for fac>i centered cubic nickel which was computed 

A a) 0 l 

by Sno\^ and Waber assuming a Dirac-Slaccr SCI* atomic calculation with 3d 4s 

configuration as input. The other parameter was Slater's exchange coefficient 

45 

ct^, which was used as unity in accord with Slater's original Xor derivation . 

ORiGiNAL P^E IS 
12 OF POOR QUALITY 



Ki^ a louct value appvoachlng tvo»Vhlrd« boon uaod, the d-b^nds vould have oppeared 

at more poiitive vnluea on this energy icaU. It la important to point out that 

the vacuum level is a natural one £or comparing a variety o£ calculations* That 

it» the potentials for either a cluster calculation or a band calculation 

depend on luperimposition SGF atomic potentials, For these, the zero corresponds 

to an electron removed to infinity* Follov/ing superposition^ this level is not 

modified* Ihe Xa procedure of sumning the individual SCF charge densities 

before taking the cube root doos not alter this potential scale since the net 

441 

charge density vanishes at *l*» for a semi-infinite solid. Snow and Waber ' 
based all their band calculations on this scale rathor than comparing the E(k) 
values with the muffin-tin potential (taken as zero.) 

* # 

In the right-hand side of Fig. Sa, the variation of the molecular orbitals 
in a cluster of oxygen approaching a cluster of nine nickel atoms* are 
8h<t».rn Th<i nicn <j to thc distance from the oxygen *u.om to cnc 

of the plane passing through nickel nuclei. Both the doubly degenerate 
E and non-degenerate A. level falls rapidly as the potentials from the nickel 
atoms increase the potential seen at the oxygen nucleus. This leads to more 

i 

negative, eigenvalues. Transfer of an electron from the nickel towards the oxygen 
and bond formation would be anticipated. So far, the perturbation of the nickel 
atoms in the cluster caused by the o:<ygcn have not been extracted from the 
calculated data. 

■ . 

Comparing the levels with only five nickel atoms in Fig. I with the 
right-hand side of Fig. 5a, the possibility of transfer is even stronger since 

*Avcril?, Adachi 4 Ellis^^^ assumed free atom configuration of 3d®4s^ for 
nickel in these clusters. 
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OF POOR QUALITY 

th«sc Ni-only valued lie at more positive velucs. The value of a that Averill 
et al used was 0.70843 and tho separation between Ni atoms was 4.709 Bohr. 

The work function for pure nickel varies between 4.07 and 5,1 eV according to 

47 \ 

Hagstrom ' but one has to take into account the transition state effects associ- 
ated with promotion of an electron into the continuum and creation of a hole in 
the d-bands. It v;ould seem, not worthwhile in the preliminary stage before achiev- 
ing self-consistency in these SSO calculations to argue about exact locations. 

I *> 

If indeed tho Fcrr.‘.l level in nickel metal is more positive, the importance of 
this work to chemisorption remains true a fortinrt . 

> 

A similar situation obtained Fig, 5b--thc right-hand side pertains to the 
shift in the molecular orbitals of carbon monoxide as it approaches a 
cluster of 5 nickel atoms, This graph has been shifted upward by 0.2 Ry to show 
the lov;er lying states in COr however, the vacuum level is no Ipncer directly 
indicated. As above^ the presence of the nickel atoms gives rise to an increas- 
ingly negative potential experienced more by the C atom than the Q atom. Levels 

shift suf ficienily to facilitate electron transfer and bond formation. 

1 

These results appear to be the first time that it has been shown that the 
metal is directly influencing either atomic or molecular levels as species 
approach a surface. Gadzuk and Schrieffor discuss such transfer processes and 
resonance broadening in terms of the free atom levels (plus ionization potentials 
and electron affinities) in the approaching species. The present results show 
that the situation is more complicated than indicated in previous treatments. 

I 

The use of the density of states is pertinent for chemisorption since 
relatively large specimens and particles have been used experimentally. The 
work on small clusters is pertinent to general problems in catalysis since 
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•mall aggregate's o£ dimensions of 10 to 5U X arc dispersed on a support such as 

amorphous silica or alumina, The"Fcrmi level*,’ below which the states arc not 

occupied in the nickel clusters, is for and for T^, namely, in t\ic 

range of 0,3 to 0,4 Ry. Note that this level is reasonably consistent with the 

. 48) 

experimental value of 5cV for the work function. The work by Johnson and Messmer 
on a single nickel cation surrounded by six oxygen anion o” is related to the 
results in Fig, 5a. They obtain a similar ’'Fernni level." 


The results in Fig. should bo compared with the results which iicssmer, 
Yong and Johnson'^^^ are currently obtaining for CO and ethylene on small metallic 
clustery. Doubtless similar trends will be found. 


* 47) 

Note* that the surface field analysis by Rao and Waber ^ has indicated some 

of the trends found in the present wore detailed molecular analysis. In Fig, 

3a of their paper, for a as for fhc oxygen ecem., tr.t 

and Py levels were degenerate and are energetically above p^. The relative 

shifts of the two levels was small (ca, 0,07, ev) but their positions "flipped" 

as the atomiSi. approached nearer the plane of nickel atoms, than n/a of 0,6, In 

reference to their Fig. 10, the barycenter these two levels was strongly 

shifted and depends on the evaluation of the Yq spherical harmonic which is not 

usually tabulated in crystal field calculations. The present Fig. 5a clearly 

shows the predicted shifts. The doubly degenerate E levels lie above when 

the oxygen atom is sufficiently far away but occurs below their barycenter ’when 

the oxygen separation is less than three Bohr units. The shift in the barycenter 

is an order of magnitude larger than the separation. 
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Fig. I* The comparison of tho molecular oigcvwaluoa for eXuscers of five 

nickel atoms ariMnsoU with and symmetry. Result of the first 
Iteration with the MSXa or "Muscatel" computer program. 

■0 

Fig. 2. Similar Comparison of the Hybrid States derived from more tightly bound 
electrons, 

> a 

Fig. 3a, Comparison of the results Cor states associated with 3d electrons 

j obtained with different computing techniques. The second column is for 
D^T-I-SSO calculations with a spherically averaged interatomic potential 
so that they can be compared with the Muscatel values. The third 
column illustratod the effect of relaxing this smoothing procedure. 

Both C\’:i*SSO calculations were done with H»50D, 

Fig. 3b. The gradual shifts in the levels for the cluster of five nickel 
atoms with increased number N of sampling point. The last column 
Illustrates the effect of Increasing the interatomic distance from 4.709 
to 6.5 Bohr units. 

Fig. 4. Illustration of the shift in individual energies of molecular states 
with successive iterations towards self-consistency. 

Fig. ’5a. Composite graph showing on the loft, the density of states cun/e derived 
by Waber and assunsing face center cubic nickel and a Slater's 

exchange coefficient of a »1. The Fermi level appears to be lower 
than the observed work function of about 0.35 Ry. On the right, the 
shift in the molecular eigenvalues cC an oxygen atom ns it approaches a 
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cLuscer of 9 nlckol atoms. 
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0 

A similar comparison showing the downward drift of molecular states 
In carbon nonoKido as it approaches a cluster of 5 niclccl atoms, 

t 

This has boon shifted upward by 0.2 Ry, compared Co Fig, Sa^ to show 
additional CO states. 
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